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' Abstract 

We have investigated the model of dark energy interacting with dark matter by choosing inhomo- 
geneous equations of state for the dark energy and a non-linear interaction term for the underlying 
, interaction. The equations of state have dependencies either on the energy densities, the redshift, the 

£>■ ' Hubble parameter or the bulk viscosity. We have considered these possibilities and have derived the 

, effective equations of state for the dark energy in each case. 
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One of the outstanding developments in astrophysics in the past decade is the discovery that the expansion 
of the universe is accelerated, supposedly driven by some exotic vacuum energy [H [2| El SI [5] . Surprisingly, 
the energy density of the vacuum energy is two-third of the critical density ( £l\ ~ 0.7) apart from dark 
matter (f2 m ~ 0.3). The astrophysical data suggest that this change in the expansion history of the 
universe is marginally recent (z ~ 0.7) compared with the age of the universe. The nature and composition 
of dark energy is still unresolved, but by using thermodynamical considerations, it is conjectured that 
the constituents of dark energy may be massless particles (bosons or fermions) whose collective behavior 
resembles a kind of radiation fluid with negative pressure. Moreover, the temperature of the universe 
filled with dark energy will increase as the universe expands [6]. The earliest proposal to explain the 
recent accelerated expansion was the cosmological constant A represented by the equation of state (EoS) 
p = —p (or u> = —1) having a negative pressure. In order to comply with the data, the cosmological 
constant has to be fine tuned up to 56 to 120 orders of magnitude [7], which requires extreme fine tuning 
of several cosmological parameters. It also posed the coincidence problem in cosmology (the question of 
explaining why the vacuum energy came to dominate the universe very recently) [8] . This latter problem 
is addressed through the notion of a tracker field Q, in which the tracker field rolls down a potential 
V(Q) according to an attractor-like solution to the equations of motion [9]. But here the field has 
difficulties reaching co < —0.7, while current observations favor uj < —0.78 with 95% confidence level [TO] , 
It is shown that a quintessence scalar field coupled with either a dissipative matter field, a Chaplygin 
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gas (CG) or a tachyonic fluid solves the coincidence problem [TT]. These problems are alternatively 
discussed using anthropic principles as well [12j . Several other models have been proposed to explain 
the cosmic accelerated expansion by introducing decaying vacuum energy [131 1 14] , a cardassian term in 
the Friedmann-Robertson- Walker (FRW) equations [15], a generalized Chaplygin gas (GCG) [16] and a 
phantom energy (u < —1) arising from the violation of energy conditions [17], [TBI US] ■ Another possibility 
is the 'geometric dark energy' based on the Ricci scalar R represented by = R/12H 2 , where H is 
the Hubble parameter [10]. Notice that 3? > 1/2 represents accelerated expansion, and 5i > 1 gives a 
super-accelerated expansion of the universe, whereas presently -R = 1/2. 

Models based on dark energy interacting with dark matter have been widely investigated [20] EH [22] 
[23] Ell Ell [Ml E3 Ell E3 EQl EIJ . These models yield stable scaling solution of the FRW equations at late 
times of the evolving universe. Moreover, the interacting CG allows the universe to cross the phantom 
divide (the transition from ui > — 1 to u> < —1), which is not permissible in pure CG models. In fact it is 
pointed out that a phantom divide (or crossing) is possible only if the cosmic fluids have some interaction 
[32j . It is possible that this interaction can arise from the time variation of the mass of dark matter 
particles [33] . It is shown that the cosmic coincidence problem is fairly alleviated in the interacting 
CG models [34j . This result has been endorsed with interacting dark energy in [35 . There is a report 
that this interaction is physically observed in the Abell cluster A586, which in fact supports the GCG 
cosmological model and apparently rules out the ACDM model [36]. However, a different investigation 
of the observational H(z) data rules out the occurrence of any such interaction and favors the possibility 
of either more exotic couplings or no interaction at all [37j . In this context, we have investigated the 
interaction of the dark energy with dark matter by using a more general interaction term. We have 
focused on the inhomogeneous EoS for dark energy as these are phenomenologically relevant. 

The outline of the paper is as follows. In the next section, we present a general interacting model for 
our dynamical system. In the third section, we derive the effective EoS for the interacting dark energy 
by employing different inhomogeneous EoS, having dependencies on various cosmological parameters. 
Finally, we present our conclusion. 

2 The interacting model 

We assume the background to be a spatially homogeneous and isotropic FRW spacetime, given by 

dr 2 



ds z = -dt z + a z (t) 



1 — kr 2 



+ r z (dQ z +sin z 6d<j) z ) 



(1) 



filled with the two component fluid namely dark energy and dark matter. Here a(t) is the scale factor and 
k = —1,0, 1 represents the spatially hyperbolic, flat or closed universe, respectively. The corresponding 
Einstein field equation is 



a \ k k 



H = - =oP-~, (2) 
\a J 3 a z 

where n = 8nG and p = p\ + p m . Moreover, the energy conservation for our gravitational system is given 

by 

p + '3H(p + p) = 0, (3) 

where p = p A = ojapa and p m = or u m = 0. We assume a special form of the interaction Q = Tp\ 
between dark energy and dark matter, where T is the decay rate. Then Eq. (3) can be subdivided into 
two parts, corresponding to p\ and p rn as follows: 

p A + 3H(l + cj A )p A = -Q, (4) 

p m + ZHPm = Q, (5) 
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respectively. Eqs. (4) and (5) show that the energy conservation for dark energy and matter would not 
hold independently if there is interaction between them but would hold globally for the whole interacting 
system as is manifest in Eq. (3). We further define the density ratio r m , by a form of scaling relation, by 
r m = pm/pK- To study how this density ratio evolves with time, we differentiate r m with respect to t: 

dTj n p m 

at pa 

Using Eqs. (4) and (5) in (6), we get 

l + r m r 



Pm 
Pm 



PA 
PA 



(6) 



= 3Hr n 



r m 3H 

Furthermore, we define an effective EoS for dark energy and matter by [38) 



eff 

J A = 



r 

3#' 



eff 



i r 

m 3H 



(7) 



(8) 



which also involve the contribution from the interaction between matter and dark energy. Using Eq. (8) 
in Eqs. (4) and (5), we get 

(9) 
(10) 



PA 
Pm 



+ 3i?(l + c4 /j > A = 



3H(l + u, 



■ff\ 



= 0. 



From the standard FRW model, the density parameters corresponding to matter and dark energy are 
defined by 

a m = f^-, = (ii) 



The above parameters are related by Q r , 



Pm n Pa 

n = , "A = • 

Per Per 

£Ia = E Using the definition of r m , we can write 

_ 1 ^A 



r^A ^a ^ 

The value of r m decreases monotonically with expansion and varies very slowly at the present era. 
Contrary to the noninteracting decreases slower when there is an interaction [3^. Using the 

definition of r m in Eq. (2), we get 

ff 2 = |(l + r m )pA- (13) 



Next, we choose the following generalized interaction term 

Q = ZHc[ 1Pm + (3 PA + 5{ PmP A) 1/2 ] 
with the corresponding decay rate 

r = 3Hc(f3 + 7 r ro + SJr n 



(14) 



(15) 



which follows from Q = Tp\ where n, /3, 7 and S are constant parameters. The above-mentioned c is the 
coupling constant. Notice that c > yields conversion of dark energy into dark matter, and vice versa if 
c < 0. Note that for j3 = 7 = n = 1 and 6 = 0, Eq. (14) reduces to the usual linear interaction term [40] . 
Making use of Eq. (14) in (4), the EoS parameter becomes 



uja = -1 - 



PA 



3Hp A 3H 

Using Eq. (16) in (8), the effective EoS of the dark energy is given by 



eff 



-1 



P A + J_Q _ n-l) 

3Hp A + 3H [L Pa >■ 



(16) 



(17) 



In the forthcoming discussion, we shall determine the effective EoS for dark energy corresponding to 
various equations of state. 
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3 Inhomogeneous equations of state for dark energy 



The general EoS relating the pressure density p and energy density p is given by 

F( Pa ,Pa)-0 (18) 

We shall also consider equations of state depending on either the redshift z, the scale factor a(t) or the 
bulk viscosity £. Note that we are not considering an EoS explicitly containing the time t, as we always 
have the opportunity to use a(t) = et~ x , with e and A constant parameters. 

3.1 Generalized cosmic Chaplygin gas 

We now take the EoS of the generalized cosmic Chaplygin gas given by [43] 

p^ = -pl a [C+( P 1 + a -C)-°l (19) 

where C = — 1 with a > 1 and A constant parameters and — I < a < 0, where I > 1. This EoS 
reduces to that of the generalized Chaplygin gas if furthermore a = and to the Chaplygin gas if further 
a = 1. The motivation to use this EoS is to construct the cosmological models that are stable and free 
from nonphysical behaviors even when the vacuum fluid behaves like a phantom energy [44) . 

Using the energy conservation principle, the density evolution is 

p A = [C+(l + Cia~ 3(1+Q)(1+ff) )TT?]T^, (20) 
where C\ is the constant of integration. We define 

Aj = {p\ +a - C) 1+a - 1 = Cia- 3(1+Q)(1+C7) . (21) 
Making use of Eqs. (20) and (21) in (17), the effective EoS for dark energy becomes 



eff_ a^i + aq— r 



A 



C + (1 + Ai)- 
where p\ is determined by Eq. (20). 

3.2 New generalized Chaplygin gas 

Zhang et al. [33] suggested another general form of Chaplygin gas, called the new generalized Chaplygin 
gas given by 

PA = Z4^), A(a) = -to A Aa-^+^ 1+a \ (23) 
Pa 

Here a is a constant parameter. This model is dual to the interacting XCDM model, where the X part 
corresponds to quintessence or X-matter (uja < —1/3). 

In this model, the energy density evolves as 

p A = lAa-^+^l^ + C 2 a-^ 1+a ^ 1+a \ (24) 
where C2 is a constant of integration. Thus using Eq. (24) in (17) the effective EoS is given by 

^^-i + ^^ + ^d-pr 1 ), (25) 

where 

A 2 = A + C 2 a 3uJAil+a \ (26) 

and pa is determined from Eq. (24). 



4 



3.3 Generalizing the generalized Chaplygin gas 



Sen and Scherrer [45] suggested an EoS for the generalized Chaplygin gas by assuming the constant 
parameter a to be free, where we have 

" A = ~A s + (l-A s )(f-)-^+a)> (27) 

where 

^ = -4^- (28) 

Using Eq. (27) one can have various cosmological scenarios: for < A s < 1 and a > — 1 we have the 
standard generalized Chaplygin gas model giving dark matter-dark energy unification. For A s > 1 and 
a > — 1, it gives the early phantom generalized Chaplygin gas, i.e., it behaves as phantom energy at early 
times and behaves like the cosmological constant lj a = — 1 at late times. For < A s < 1 and a < — 1 it 
represents the transient generalized Chaplygin gas, in which case Eq. (27) gives de Sitter regime at early 
times and a matter dominated regime at later times. 

The density evolution is given by 



Pa = Pa 
The corresponding effective EoS is 



l/l + Q 

A. 5 + (i-A,)(-r 3(1+Q) 

a„ 



with pa determined by Eq. (29). 

3.4 Interacting scale factor dependent dark energy 

We here take the EoS [46] 

Pa — — Pa(1 + Aa a ). (31) 

The corresponding density evolution is 

p A = C a expf— \ (32) 



(29) 



V a 

with C3 is constant of integration. The effective EoS is given by 

w e// = _ j4aQ + ^ (1 _ p „-l ) (33) 

3.5 Interacting Hubble parameter dependent dark energy 

An interesting EoS depending on the Hubble parameter H is given by [47) 

PA = -PA + .f(pA) + G{H). (34) 

The corresponding FRW equation is 

PA = -3H[f( PA ) + G(H)}. (35) 

Let us choose the following EoS: 

f( PA ) + G(H) = -Apl - BH 2 \ (36) 
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where e is a constant. Using Eq. (13) in (36), we get 

f( PA ) + G(H) = -Apt - B'pl, (37) 

where 

B' = B 

.3 



B 1 = B 

Using Eqs. (35) and (36) in (17), we get 



?(l + r m ) • (38) 



ojf s = -1 - (Apt' 1 + B'pX 1 ) + ^-(1 - Pr 1 )- (39) 
Here p\ is determined from Eq. (13). 

3.6 Interacting redshift dependent dark energy 

We here assume that the dark energy evolves with the redshift parameter z. Hence we take the following 
linear EoS 0S]: 

u(z) = u + uiz, (40) 

where lo q and io\ are constants. This EoS was used to analyze the cosmic microwave background and 
the matter power spectrum |49j . Eq. (40) effectively works for z < 1, while lo(z) = lj q + uj\ holds up till 
z 1. The thermodynamical properties of dark energy have been investigated using Eq. (40), and it is 
deduced that the apparent horizon of the universe may be the boundary of thermodynamical equilibrium 
for the universe like the event horizon for a black hole [50]. We are interested in the evolution of dark 
energy (i.e. Eq. (40)) in our generalized interacting model. The energy conservation principle gives 

PA = C 3 a-^°- i)exp(3 Wl ^), (41) 

where we have used z = (a /a) — 1 and C3 is a constant of integration. Notice that for u>i — 0, Eq. (41) 
gives the evolution of the usual dark energy. 

Using Eq. (41) in (17), we get 

cul" =oj + oj 1 z + ^(1-pI- 1 ). (42) 

h0.08 



The astrophysical data support cosmologies with ui = —1.25 ± 0.09 and u>i — 1.97_q 01 E3- Also using 



< c < 1, we see that the right hand side of Eq. (42) becomes negative, i.e. w e A < 0, thereby supporting 



the existence of phantom energy. 
We now take another EoS PUIIST] : 

u>(z) = uj + uji ^1 - = lo + Wx ( T~^ J ■ (43) 

It avoids the divergent behavior as opposed to Eq. (40) and hence is used to parameterize the astrophysical 
data to higher redshift, up till z « z rec . The two constants appearing in Eq. (43) are constrained: 
-1 < lo < -0.434 and -0.564 < u>i < 0.498 52J. Its density evolution is given by 



PA 



where C4 is a constant of integration. Thus the effective EoS is 



= C 4 a" 3(1+ ^ +Wl) exp f 3wi— J , (44) 



u B /t = u + W i ( 1 - 3^) + ^(1 - p\- x ), (45) 



with pa is determined from Eq. (44). 
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3.7 Interacting viscous dark energy 

The Eckart theory 53J effectively deals with fluids having non-zero viscosities. The term viscosity arises 
from fluid dynamics, which has two major parts, namely the bulk viscosity £ and the shear viscosity 
77. In viscous cosmology, shear viscosities arise in relation to space anisotropy, while the bulk viscosity 
accounts for the space isotropy [SHEB, 56J . The bulk viscous fluid is represented by the Eckart expression 
II = — £(p A )utt, where is the four velocity of the viscous fluid. The bulk viscosity is generally taken to 
be positive to ensure positive entropy production in conformity with the second law of thermodynamics 
[57] , Its scaling may be represented by £ ~ pT^ , where C is a constant parameter. 

The effective pressure containing the isotropic pressure and the viscous stress is given by 

Peff — PA + II, (46) 

where p\ — x/ Pa with X > an d II = — 3-ff£ in the FRW model. Thus the energy conservation principle 
for the bulk viscous fluid becomes 



PA + 3H{p A + PA -3HZ) = 0. 



(47) 



Using Eq. (13) in (47) we have 



Pa + 3H(p A +p\- S.Yy/PA) = 0, 



where T = y%3K(I + r m ). Solving Eq. (48) we get 



a(t) 



exp 



P A +a ^^PA+X 



-1/3 



(48) 



(49) 



1/2 

with C5 is a constant of integration. This equation can be solved exactly by choosing £ = vp A with v 
a constant; thus, we have 



PA 



(C 5 a 



-3(l-uT)(l+aO 



X 



1 - vT 



Defining 



A3 = (<7 5a )-3(l-«T)(l +01 ) = p l+a {1 _ vT) + ^ 

Using Eqs. (50) and (51) in (17), we get the effective EoS of the interacting viscous dark energy: 



u> 



l ff - -1 + (1 - vT) 



X) 3H 



(50) 
(51) 

(52) 



4 Conclusion 



In this work, we have determined various effective equations of state for the dark energy having non- 
zero interaction with the matter in the universe. The dark energy can have dependencies on various 
cosmological parameters like the Hubble parameter, the redshift, the scale factor, the energy densities or 
the bulk viscosity. We have considered all such possibilities in our interacting model. 
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